Global scale-invariant dissipation in collisionless plasma turbulence 
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A higher-order multiscale analysis of the dissipation range of collisionless plasma turbulence is 
presented using in-situ high-frequency magnetic field measurements from the Cluster spacecraft 
in a stationary interval of fast ambient solar wind. The observations, spanning five decades in 
temporal scales, show a crossover from multifractal intermittent turbulence in the inertial range to 
non-Gaussian monoscaling in the dissipation range. This presents a strong observational constraint 
on theories of dissipation mechanisms in turbulent collisionless plasmas. 
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The solar wind provides an ideal laboratory for the 
study of plasma turbulence [T]. In-situ spacecraft obser- 
vations suggest well-developed turbulence at 1 AU with 
a magnetic Reynolds number ~ O (lO 5 ) [2] [3]. These 
show an inertial range of Alfvenic turbulence on mag- 
netohydrodynamic (MHD) scales which is an anisotropic 
and possibly compressible energy cascade |4j |5j [6] with 
intermittent magnetic field fluctuations described by sta- 
tistical multifractals and a power spectral density (PSD) 
with a scaling exponent close to —5/3 [T]. An outstand- 
ing problem is how, in the absence of collisional viscos- 
ity in the solar wind, this inertial range of MHD turbu- 
lence terminates at smaller scales where one anticipates a 
cross-over to dissipative and/or dispersive processes via 
wave-particle resonances. Understanding the nature of 
the dissipation processes may also inform open questions 
such as how the solar wind and solar coronal plasmas are 
heated [HI [Sl- 
it has long been known [TUJ, [TT] that in collision- 
less plasmas there is a transition in the PSD at high 
wavenumber k from MHD to kinetic physics at approx- 
imately the ion gyroscale. High resolution in-situ mag- 
netic field observations reveal that at these scales the 
turbulent solar wind shows a transition from a ~ —5/3 
power law in the inertial range to a steeper power-law at 
higher k with spectral exponents in the range (—4, —2) 
[T21 MS] . However, the relevant physical mechanism is 
much debated; having implications for phenomena as di- 
verse as magnetic reconnection [141 115) . neutron stars 
and accretion disks [T5]. Theories which have been pro- 
posed range from nonlinear turbulent-like cascade pro- 
cesses [T7J HHJ US] to weak turbulence theories with wave 
dispersion and resonant plasma interactions [20]. As well 
as studies of in-situ spacecraft measurements in the solar 
wind [2T] , foreshock [22] and magnetosheath [23] [23] re- 
gions, these theories are explored using simulations rang- 
ing from Hall-MHD [25], electron-MHD [HI [25], gyroki- 
netics [27 , particle-in-cell simulations of whistler turbu- 
lence [28] and Vlasov- hybrid simulations [29] . 

Both neutral fluid and MHD turbulence share a 'clas- 
sic' statistical signature - namely an intermittent mul- 
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Figure 1: PSD plots of the components of the magnetic field 
from both FGM (at frequencies lower than 1 Hz) and STAFF- 
SC (at frequencies above 1Hz) instruments. The PSD values 
for B x and B y have been shifted up for clarity. The 95% confi- 
dence intervals for all these plots are at ±0.03 log^nT 2 H z" 1 . 
The shaded region denotes the frequency range of STAFF-SC 
to be studied; the upper limit corresponds to the search-coil 
signal-to-noise ratio (SNR) of lOdB; the lower limit to the 
SNR of 20dB [30]. At lower frequencies the STAFF-SC re- 
sponse is attenuated by the instrument callibration filter. 



tifractal scaling seen in the higher-order statistics. In 
this letter we test the statistical properties of the dissi- 
pation range and find in contrast monoscaling behaviour 
i.e. a global scale-invariance. This provides a strong 
discriminator for the physics and phenomenology of the 
dissipation range in collisionless plasmas. 

We present a detailed analysis of an interval of quiet, 
stationary solar wind observed in-situ by the Cluster 
spacecraft quartet during an hour interval 00:10 - 01:10 
UT on January 30, 2007. By combined analysis using 
high-frequency measurements of magnetic field fluctua- 
tions from the search-coil (STAFF-SC) [3T] and flux-gate 
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Figure 2: Summary plot illustrating the quiet nature of the 
solar wind interval under study. Top panel: E-field spectro- 
gram from WHISPER on Cluster spacecraft 4 showing steady 
plasma emissions; center panel: B-field components; bottom 
panel: ion velocities from HI A on Cluster spacecraft 1. 



magnetometers (FGM) [32) , we simultaneously probe 
both the incrtial range and dissipation range up to fre- 
quencies of 80 Hz. This will allow us to establish the 
scaling properties in the dissipation range over about two 
decades in frequency. Figure [T] shows the overlaid PSD 
from FGM and STAFF-SC for our interval for the three 
magnetic field components B Xl B y and B z in geocentric 
ecliptic (GSE) coordinates - both the inertial and dissi- 
pation ranges can be clearly seen in all three components. 
We indicate on the plot the frequency corresponding to 
the Doppler-shifted proton gyroradius f pg [TH] which is 
close to the cross-over frequency break /f, r between the in- 
ertial and dissipation range. A summary of the solar wind 
interval used in this study is shown in fig. [2] We plot the 
electric field spectrogram in a frequency range 2-40 kHz 
from WHISPER [33 , the magnetic field timeseries from 
FGM and the ion velocity from CIS/HIA [34 during the 
interval of interest. The constant ion velocity and plasma 
density (evidenced by the constant electron plasma fre- 
quency on the E-field spectrogram) indicate a stationary 
pristine interval of ambient solar wind disconnected from 
the Earth's bowshock, away from ion and electron fore- 
shock regions and imbedded in a steady fast solar wind 
stream > 650 kms^ 1 i.e. quintessential plasma turbu- 
lence away from any external a priori physical processes 
and dynamics. Due to stationarity of plasma parameters 
and the B-field magnitude it is sufficient to quote single 
values for the other relevant parameters: ion tempera- 
ture Ti ~ 1.2 MK, Alfven speed Va ^ 50 ferns -1 , ion 
plasma beta j3 ~ 2, and plasma density n e ~ 3.8 cm -3 
(from the plasma frequency measured by WHISPER). 
Importantly for this study, we have chosen an inter- 



val where both STAFF-SC and FGM were operating 
in burst-mode so that one can access the largest range 
of scales available and probe further into smaller scales 
where the dissipation range can be studied. For this in- 
terval STAFF-SC provides AC waveform data in a fre- 
quency band between ~ 0.1 Hz and 180 Hz (low-pass fil- 
ter); our study is restricted to frequencies lower than 80 
Hz to maintain good SNR. FGM provides DC waveform 
data with largest frequency at ~ 33 Hz (Nyquist cutoff). 
This ensures a large sample (~ 1.6 x 10 6 for STAFF-SC 
and^ 2.5 x 10 5 for FGM) and thus well-resolved statis- 
tics. The Welch PSD method used in fig. [I] employs 
50% overlapped windows each containing a sample of 2 14 
points for FGM and 2 16 points for STAFF-SC; this pro- 
vides a very large range of frequencies to study as well 
as reducing the error due to noise on the PSD measure- 
ments. The large sample size coupled with the intervals 
stationarity ensures that we have a good control of errors 
which can arise due to finite sample size |35j . All the re- 
sults presented in this paper are from Cluster spacecraft 
4; the B z component provides the largest overlap in fre- 
quencies between both FGM and STAFF-SC due to it 
being out of the spin plane of the spacecraft and thus is 
cleanest with respect to spin tone contamination. As we 
will be studying single point Eulerian measurements in 
very fast solar wind, streaming past the spacecraft pro- 
viding a single time-series, we will assume the validity of 
Taylor's frozen-in-turbulence hypothesis [23 which uses 
time as a proxy for space - although our PSD and other 
statistics will always be presented in the frequency and 
time domains. 

The PSD provides one statistic to probe the scale de- 
pendent behaviour of the turbulent fluctuations and is 
equivalent to studying the autocorrelation - a second or- 
der statistic. To test for multifractal scaling of the fluc- 
tuations we now turn to higher-order statistics. We focus 
on the statistics of magnetic field increments defined as 
SBi(t,r) = Bi(t + t) — Bi(t) for each vector component 
i and time lag or scale r; in particular the focus is on 
the absolute moments of these increments, also known as 
structure functions 

1 N 

5 rw = ^Ei^> T )r . (!) 

where m is the moment order. We have verified the 
statistical stationarity of the interval being studied and 
thus can form an ensemble average by taking a time av- 
erage (assuming ergodicity) along the signal of sample 
size N. Importantly, the higher-order structure functions 
progressively capture the more intermittent, larger fluc- 
tuations. As we are studying the magnetic field incre- 
ments, these large fluctuations represent the spatial gra- 
dients which are responsible for dissipating energy from 
the magnetic fields. We will focus on the scaling be- 
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Figure 3: Structure functions of orders: 1-D, 2-V, 3-*, 4-o 
and 5-0- Open shapes correspond to FGM measurements 
and filled shapes refer to STAFF-SC. The curves have been 
shifted along the vertical axis to allow a comparison of the 
gradients. The shaded area indicates the scales where both 
FGM and STAFF overlap. Linear fits for the inertial and 
dissipation ranges are also shown. 

haviour of the structure functions with scale r such that 

S?(t) ex r<M , (2) 

where linear dependence of the scaling exponent C( TO ) = 
Hm implies monoscaling with a single exponent H. In 
theories of turbulence non- linear £( m ) behaviour is as- 
sociated with the intensity of energy dissipation being 
distributed on a spatial multifractal pQ. 

The structure functions for the data interval studied 
here are shown in fig. [3] On this log — log plot the gra- 
dients as shown give estimates of the scaling exponents 
C(m). The inertial and dissipation ranges of scaling are 
well-defined with a sharp transition at the break point at 
~ 3 seconds in agreement with the PSD in fig. [I] the dis- 
sipation range extends over nearly two orders of magni- 
tude. Importantly, there is excellent agreement between 
STAFF-SC and FGM in the dissipation range where they 
overlap for almost a decade, indicated by the shaded re- 
gion on the plot. We plot C( m ) vs - m f° r the dissipation 
range in the main panel of fig. |4j and for the inertial 
range in the inset. Surprisingly, the dissipation range is 
monoscaling i.e. globally scale- invariant; in contrast to 
the inertial range which is multifractal, characteristic of 
fully developed turbulence with ((2) ~ 2/3. The sin- 
gle scaling parameter for the dissipation range for B z is 
H = 0.89 ± 0.02 for STAFF-SC and H = 0.84 ± 0.05 for 
FGM. To test the robustness of this result we have re- 
peated this analysis for another ambient fast solar wind 
interval (12:10 - 14:00 UT January 20, 2007) and ob- 
tained the same global scale-invariance. In both cases we 
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Figure 4: Main plot: Scaling exponents C, with order m; a lin- 
ear relationship on this plot indicates monoscaling behaviour. 
C(m) obtained from both FGM and STAFF-SC are shown 
for B z ; these show close correspondence. STAFF-SC B x , B y 
components are also shown and indicate isotropic scaling. In- 
set: £(m) Vs. m for the inertial range using FGM B z \ this is 
concave, consistent with the multifractal nature of the inertial 
range. 

find that all three field components are monoscaling. We 
can see that for the particular solar wind interval shown 
in fig. [4]the exponents H for B x and B y are close to that 
of B z , suggesting that the small scale features of this tur- 
bulent interval of the solar wind are also isotropic. For 
STAFF-SC data from the second interval, however, we 
find H = 0.9 ± 0.02 for B x and B y and H = 0.8 ± 0.05 
for B Zl suggesting an anisotropy that may depend on 
local plasma parameters. 

Monoscaling of the structure functions implies that 
the probability density function (PDF) of the increments 
P(SBi,r) at a particular scale r should collapse onto 
a unique scaling function V s via the following rescaling 
operation [36 

V s (SB 1 t- h )=t h P(SB 1 ,t) . (3) 

This collapse of the data to a single scaling function is 
tested in fig. [5] for B z , where we have used the same val- 
ues of r and the H value obtained above. We can see that 
there is an excellent collapse onto a single curve. A fit- 
ted Gaussian illustrates the highly non-Gaussian nature 
of the tails of this PDF. 

In conclusion, our results suggest that global scale in- 
variance in small-scale magnetic fluctuations is a robust 
feature of the dissipation range of collisionlcss plasma 
turbulence in the fast ambient solar wind. This is a sur- 
prising result as it is distinct from the multifractal scaling 
that is characteristic of both neutral fluid and MHD tur- 
bulent cascades in the inertial range. Successful theoreti- 
cal understanding of the dissipation range should include 
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Figure 5: Main plot: PDFs rescaled using eq. Q 
[SB" = 5B z T~ H y Inset: PDFs at different scales r before 
rescaling; red and green curves show the smallest and largest 
values of r respectively. A Gaussian fit to the data illustrates 
the heavy-tailed non-Gaussian nature of the rescaled PDF. 



this property. Our result provides a strong discriminator 
of the relevant physics and phenomenology; for example 
the monoscaling that we find is reminiscent of that found 
at higher orders in electron-MHD simulations [37] . To de- 
termine whether this phenomenology is in fact universal, 
future studies should aim to reproduce and/or break this 
result in more dynamic environments such as at planetary 
shocks [22] > magnetosheath [24] and at sites of magnetic 
reconnection [TH [15] ; although the main difficulty here 
will be to identify sufficiently long stationary intervals. 
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